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Abstract 

Diagrammatic techniques are well-known in the calculation of dispersion interactions between 
atoms or molecules. The multipolar coupling scheme combined with Feynman ordered diagrams 
\sO significantly reduces the number of graphs compared to elementary stationary perturbation the- 

ory. We review calculations of van der Waals-Casimir-Polder forces, focusing on two atoms or 
molecules one of which is excited. In this case, calculations of the corresponding force are notori- 
ous for mathematical issues connected to the spontaneous decay of the excitation. Treating such 
Oh unstable states in a full non-equilibrium theory provides a physical interpretation of apparent 

+1j contradictions in previous results and underlines the importance of decay processes for the in- 

termolecular potential. This may have important implications on reactions in biological systems, 
where excited states may be relatively long-lived and the resonant intermolecular force may result 
in directed Brownian motion. 
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1 Tools from Non-equilibrium Field Theory in Molecular QED 
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t— I The interaction between electrically neutral yet polarizable particles, commonly named after J. D. 

van der Waals, F. London, H. Casimir, and D. Polder, is one of the fundamental problems in atomic 
and molecular physics. While the two-body potential between two particles (atoms, molecules, 
nanoparticles) in their internal ground state is unambiguous [TJ-J3J , apparently incompatible results 
have been obtained if one of the atoms is prepared in an excited state. The quantity of interest is 
the long-range part of the potential that overwhelms the familiar van der Waals interaction: some 
calculations found this part to oscillate spatially (5H7), while later work found a monotonic power 
law [8]. It is the aim of the present paper to understand these differences [9[[l0] better. We handle 
this non-equilibrium problem with a diagrammatic (Feynman graph) expansion. 

The starting point of our quantum mechanical description is the electric dipole coupling Hamil- 
tonian 

H AF = -Y, Ei{x n ) [d%r{Xn)^ 9 \x n ) + dT riXnWHxn)] ■ (1) 
n 

For simplicity, we treat the atoms as pointlike two- level objects located at the spacetime coordinates 
x n = {r n ,t n }. The fermionic operators ij) a ^(x) create atoms in state a = g,e at x. Assuming that 
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the transition dipoles are real, we work in the following with an effectively scalar electric field 
E{x) = d™Ei(x) whose space argument accounts for possible differences in the transition dipoles. 
The textbook treatment of dispersive atom-atom interactions (4 12 employs stationary fourth- 



order perturbation theory in Haf- The summation over all possible intermediate states leads to a 
rather large number of terms. In this work, we discuss a more concise treatment using the closed 
time-path contour formalism due to Schwinger, Craig, Mills, and Keldysh (see Ref. |13|). Following 



Sherkunov 10 11 14 , the level shift of a ground-state atom is extracted from its full propagator 

i9^\x,y) = (T c {S 2 (-oo, 00)^(00, -«>)^(x)^ t (y)}> conn . (2) 

The greek subscripts take the values 1 or 2, denoting the branch of the Keldysh time contour which 
extends from t = —00 to 00 (branch 1) and back (branch 2). The prescription T c {. . .} orders 
operators according to their time arguments on the contour. Finally, S is the standard scattering 
operator, the brackets (...) denote an eigenstate of the non- interacting theory , and only connected 



products contribute. The Feynman formalism 15,16 is valid for known in- and out-states, which is 



the typical scenario in equilibrium systems at zero temperature. If, however, a general initial state 
(including, in our case, unstable excited atomic states) is prepared and then left to evolve under the 
interaction, the Keldysh formalism is necessary to describe all energetically allowed processes. 

In the following, a propagator with both times on the forward branch will be used. Its series 
expansion is from Eq. Q 

ig { ^\x,y) = (T{r(x)^(y)}) 

i r 2 

- - / dhdt 2 (t c { j2 (-ir +fi HAFAti)HAFA^f(^f\y)}r nn + ■■■ ■ (3) 



where the first line is the ordinary (bare) Feynman propagator given in Eq. (16), with T being the 
standard time ordering. We will see below that processes involving only ground-state atoms can be 
entirely described in terms of Feynman propagators. 

The next section presents the calculation of self-energies for atoms and photons, and recovers 
the energy shifts and decay widths in a system of two particles. A comparison of the Feynman 
and Keldysh results suggests an interpretation of the apparent disagreement of earlier results: these 
may apply in different time regimes after preparation, separated by the time scale that characterizes 
energy transfer by an exchange of resonant photons. Finally, we discuss how resonant long-range 
potentials might influence some systems of biological relevance that are not commonly considered in 
the community present at this conference. 



2 Self-Energies with Dressed Photons 

Our perturbative analysis is built from two basic expressions: the self-energies of an atom coupled 
to the photon field and of a photon coupled to a (second) atom, respectively. Conceptually, this is 
closely related to the notion of dressed states (Tt| . The real and imaginary parts of the self-energy 

£n = A£-*~ (4) 

can be identified as the energy-shift and the inverse lifetime (decay width) of the particle. We first 
calculate the atomic self-energy at next-to-leading order in the coupling to a generic photon field. 
In a second step, the interatomic interaction is identified by "dressing" the photons with a second 
atom. We only have to consider terms where photons are connecting the two atoms, provided we 
assume the transition frequencies to include the single-particle Lamb-shift. 
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2.1 Single atom plus field: atom self-energy 

We start by considering the propagator of a ground-state atom B (transition frequency u B ) at the 
one-loop level. Using the Feynman rules from the Appendices and App. C of Ref. |16| , the evaluation 
of Eq. ^ yields 

g ( ff(x,x') = ,'o*o^%. ,- 

= J d 3 Xl J d 3 x 2 J du; e - <w t t -* / )^u(a;,r,ri)Ef 1 (a; > n,r2)^ii(w,r aj r') , (5) 
where the self energy for a pointlike atom is to lowest order 



B / JS/ ^\ B r d 

Sff (wm) = , ^ m =~i ir D u( u > rB ' r B)en(uj- m + w) . (6) 



The excited state propagator en is given in Eq. (17) and D\\ denotes the Feynman propagator for 



photons. For the self-energy Sf 1 of an excited atom, replace en by gu in the integral. This result 
can be brought into a more familiar form (see Ref. Il8|) by evaluating it on the mass shell, ui- m = e B , 
and writing the integral over positive frequencies only (recall that Du(uj) is even in u): 

Eff (e? ) = * jH ^ D n (u, r B , r s )< ( w ) . (7) 
Here we recover the Feynman-ordered polarizability 

that is evaluated explicitly in|Aj (For the excited two- level atom A*, replace co B -> -o; A .) Note that 
the propagators here have poles in the upper left and lower right quadrant of the complex frequency 
plane (Feynman prescription). They coincide, however, with the usual retarded response functions 
at positive frequencies, and only these appear in the complex self-energy ([7]). 

Substituting the bare photon propagator into S^{ e yields, according to Eq. Q, an infinite energy 
shift which is, of course, the unrenormalized Lamb shift of atom B. The imaginary part of S n is 



nonzero if the temperature T > [see Eq. (23)], describing the absorption rate of thermal photons. 
At T = 0, an imaginary part is found for the excited atom A*: this comes from a pole of af^(oj) in 
the upper right quadrant and gives the spontaneous decay rate Tq A = 2 Im .On (0^,1-^4, r a). 

2.2 Photon propagation in the presence of a second atom 

To evaluate the impact of a second atom (labeled A) on the self-energy of atom B, we take the 
diagram Q and replace the photon propagator by its next-to-leading order correction: 

A-^__j^A A*. ^A 

Cvwvo — > o^h^^v-^o or o/\/\^^5' v/x/>0 j 

depending on the state of atom A. This "modular" strategy is, of course, consistent with the fourth- 
order expansion of the atom propagator gu in Eq. ([2]). The bare photon propagator Dn(x,x') is 
given in |Bj and its correction to leading order recovers Eq. (28) of Ref. 10 

D^ )a (x',x) = -± dhdt 2 (T c { (- 1) W ' H AF)i (h)H AFtl , (t 2 )E! {x')Ei (x) }) c ° 

J Ll'=l 
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where a = e,g denotes the state of atom A. For the purpose of our analysis it is sufficient to restrict 
the field to zero temperature which simplifies the treatment in the frequency domain considerably 
(see[Bj). We find from Wick's theorem [l~3] 



D { S a {uo > 0,r',r) 



Dn(u, r, Y A )a a 1 ^{uj)D ll {oj, r A , r) 



0) 



where the atomic polarizabilities are given in section [Aj Note that this contains an extra term that 
is not in the form of Feynman propagators. There is no connected diagram involving D22, and at 
positive frequencies, D^w) = (see[B|). 

2.3 Two ground-state atoms 

If atom A is in the ground state (a = g), the contribution from the second Keldysh branch vanishes 
because af^(w > 0) = 0. We obtain a two-atom self-energy by substituting the dressed photon 
propagator D^ 9 into Sf 1; 



^11 





.4 
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A 
D 



■ [ ^a{^{u)a 9 l ^{u)Dii(uj,r A ,TB)Dii(u,YB,Y A ) . 

J0 27T 



Only Feynman-ordered quantities contribute, and indeed, this situation corresponds to the ground 
state of the non-interacting theory where no non-equilibrium formalism is needed. The integration 
lends itself to a rotation to imaginary frequencies, and one sees that £ff = AE" is a purely real 
two-body potential, equal to the well-known Casimir-Polder potential if the two atoms are embedded 
in free space. Note that the same calculation in traditional perturbation theory, e.g. in Ref. |2j, 
involves 12 diagrams with time-directed photon-lines rather than the two Feynman diagrams above. 



2.4 Ground-state atom and excited atom 

The additional term in the photon propagator ([9]) makes the Feynman and Keldysh results differ. 
We find the following self-energy for the ground-state atom B in the presence of an excited atom A* 



\ — a ii^){ D ^^^B,Y A )a\i{uj)Dii(u},T A ,Y B ) 
Jo 27T 

- D u (uj, r B , r A )al 2 (u})D2i(u}, r A , r B )} 



(10) 



+ 




B 

A 



[ID 



The first line is the result obtained in Refs. [5]- [7j. In our formalism, this part of the result contains 
only time-ordered (Feynman) quantities and is represented by the first two diagrams of Eq. ( 11 ). We 
will argue in the next section under which conditions this line gives already the full result. 

The second set of diagrams illustrates the terms that arise from the second Keldysh branch. There 
is one vertex beyond the (thick red) bars that corresponds to an interaction operator H A f^(x) on the 
second Keldysh branch. It will end leftmost after the contour ordering T c and can be interpreted as 
acting directly on the outgoing states. This modification of the out-state translates the equilibration 
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of the initial state (prepared at t = — oo and left to evolve under the interaction). The thick red bar 
thus illustrates the split between the branches of the Keldysh contour: it corresponds, apart from a 
prefactor, to cutting the diagrams and putting the loose ends on the mass-shell. The resulting two off- 
diagonal processes conserve energy (tree structure of the diagrams) and can be read as spontaneous 
emission from atom B after a resonant photon exchange (left diagram) and as emission from atom 
B stimulated by a photon from atom A (right diagram). 

To provide a more transparent physical interpretation, we bring the above result into a form 
involving retarded and advanced response functions. These are related to the Keldysh propagators 
by an = dR+ai2 and D21 = D\\ — D A (see Refs. (l3fl9] ). Remembering that for positive frequencies 
and T = 0, we have Du(w) = D r {oj) and af^uj) = a R (iv), we get 



^11 



J ^a 9 ^(uj){D R (uj,r B ,r A )a e R (uj)D R (uj,r A ,r B ) 

+ Dr(w, r B , TA)a.l 2 (u)D R {uj, r A , r B ) 
- D r (uj, r B , r A )al2(uj)D R (u), r A , r B ) 
+ D R (uj,r B ,r A )al2(uj)D A (u,r A ,r B ) } . 



(12) 



The two lines indicated by the asterisk (*) cancel out exactly. The first line gives an integrand 
regular in the upper right quadrant and can be rotated onto the imaginary axis; it does not yield 
any oscillating contributions. Within the two- level approximation, one has a e R = —a R [see Eq. (JsJ) ] , 
so that the first term is equal, up to a global sign, to the ground-ground self-energy Sff. The last 



line gives a purely resonant contribution because af^(cj) = 2ni5{uj — uj a ) according to Eq. (21). It 
involves the modulus squared of the (retarded) photon Green function because D A {uS) = D* r {oj) at 
real frequencies. We thus rewrite Eqn.(fl~2l) as 



^11 



_2j ll 




a 9 R B {u) A ) \D R (u A ,r B ,r A 

A A 



(13) 





A 

D 



A 
B 




A 

- B 



off (.ua) , 



where the last term can be identified as Forster resonant energy transfer (2} [12} [20]-[23] (FRET). 
Rather than distinguishing absorption and emission, the above reordering has led to a separation of 
dispersion (nonresonant or virtual photons) and FRET effects. The resonant energy exchange makes 
the ground-state atom B unstable, and we recover the rate for FRET 

Tfret = -2Im[Sf^] = 2k5((jj b - u) A )\D R (u A ,r B ,r A )\ 2 , 

in full agreement with the Golden Rule. The generalization of this rate to molecular emission and 
absorption spectra of finite width is straightforward and well-known . The FRET process comes 
along with a resonant two-body potential (the real part of £fi) that has been discussed in Refs. [8]- 
|11| . This potential does not oscillate spatially and lends itself to a simple semi-classical picture: the 
polarization energy of atom B in the time-averaged field of a dipole source located at atom A, whose 
amplitude is fixed by the transition dipole |9|. 



2.5 Transient spatial oscillations of the resonant potential 



It is instructive to evaluate the contribution that was cancelled from Eq. (12) (upper line marked 
(*)). The explicit form of af^ from Eq. (21) results in a spatially oscillating contribution to the 
self-energy 



-a 9B {uj A )Dn(u A , r B , Y A )D n (u A , r A , r B ) . 



(14) 
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Hence the first two (Feynman) diagrams in Eq. (10) give 




- A 



+ 




- A 



-Sff + A£*f 



(15) 



Since the potential £ff is monotonous, the extra term (14) is responsible for the spatially oscillating 
potential found in Refs. [5]- [7] using equilibrium (Feynman) theory. If all four diagrams are taken 
into account, one reaches a fully equilibrated state and the interaction S^? [Eq. (13)] does no longer 



contain any spatial oscillations. This is the result in Power and Thirunamachandran's calculations 
gg. ' 

We suggest the following scenario to understand the physical relevance of the two results, using 
the narrative of atom dressing 17 . The interchange of virtual (or non- resonant) photons responsible 
for the monotonous part of the potential is very fast. This short-time behaviour is similar for both 
ground-state and excited atoms. The exchange of near-resonant photons takes a longer time due to 
the small frequency differences involved, and it leads to the equilibration of the initial state. This 
happens on the scale 1/Tfret set by the Forster rate. This time appears as a lower limit for the 
validity of the equilibrated potential obtained from E^. An upper limit is set by the spontaneous 

The two time scales 



lifetime 1/Tq of the atom A* when photons appear in free space modes 11 



22 



are well separated at short distance (non-retarded range, typical for FRET experiments). 



The spatial oscillations of the potential (15) may therefore appear as a transient effect right 
after the preparation. We speculate that they remain visible if the excited state is continuously 
replenished with a sufficiently high rate (larger than Tfret)- This would require a time-dependent 
perturbative calculation similar to Ref. 12] that is beyond the scope of this paper. We note that a 
spatially oscillating dispersion interaction has been found at short times in the somewhat analogous 
situation of an excited molecule close to a surface 1241. 



3 Discussion and Conclusions 

The apparently differing predictions for the resonant part of the two-body potential between a ground- 
state and an excited atom have been a source of confusion for quite some time. We have reviewed 
this system in a non-equilibrium description (Keldysh closed time-path formalism) and identified how 
contributions that oscillate spatially disappear in such a description. The two results may correspond 
to different physical setups or to situations separated by a characteristic time scale, which we related 
to the rate of energy transfer (FRET): the oscillatory potential (Refs. [5]- J7]) holding on very short 
time scales, and eventually evolving into the non-oscillatory one of Refs. [8p |11| . 

In molecular physics, Forster broadening of atomic or molecular spectra is well known, but there 



seems to be less literature on the forces that come along with such resonant exchanges of energy 25 



Actually, the term Forster force was coined as late as 2003 by Cohen and Mukamel |22 23 . These 
forces may, however, play an important role in systems that stand in the focus of recent research, 
ranging from cold atoms [26p7] to quantum dots, NV centers, and even biomolecules such as proteins 
and DNA. In the latter context, it has been proposed that diffusion- limited reactions between an 
excited and a ground-state reactant are facilitated by a state-selective force. The force would modify 
the motion of the molecules relative to undirected Brownian diffusion which may help to understand 
unusually high reaction rates [28| . In biological systems in particular, the excited state may be 
pumped by an energy source (photoabsorption, chemical) . Then the energy flow through the system 
will determine to what extent the scenario involving equilibration after excitation actually applies. 
Such questions certainly deserve further investigation, and underline that talking about unstable 
states as in dispersion interactions between atoms or molecules requires a careful characterization of 
the excited state in question. 
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A Atom fields and propagators 

Atoms are considered distinguishable by an index n = A,B, with an internal two- level structure of 
energies e" e . The results will only depend on the Bohr frequencies uf 1 = e™ — e™ . The annihilation 
operator evolves in the interaction picture proportional to exp(— ie™t) (a = g,e). Atoms are 

considered immobile and pointlike so that the position dependence is carried by the field operator in 
the interaction Hamiltonian Q. 

The bare (leading order) propagator for a ground-state atom is given by 

fofrx) = -i(^(x')^(x))9(t' -t) = J ^ - • (16) 

For an excited atom, is obtained using e™ in place of e™. The step function 6(t' — t) arises because 
in this non-relativistic theory, there are no anti-particles. In the frequency representation, we thus 
have the Feynman rules 

5nH = 1 ^+=° °> e nH = 1 v^+ =(> "^^ ( 17 ) 

uj — e g + i u) — e e + i Ir 

We also need Keldysh propagators for the atomic polarization operator P{x n ) = ip e (x n )ip 9 ^ (x n ) + 
■4)9{x n )^(x n ) that appears in the dipole interaction The general correlation is 

a™p(x',x) = i(T c {P a (x' n )Pp(x n )}) a (18) 

for atom n in state a. 

a\l{x\x) = -{(e\^(x')^(x)\e)g^(x',x)+x' o x) (19) 

d — e— ( - - + - -) (20) 

The Feynman polarizability a^™(cj) of the excited atom has a pole in the upper right quadrant of the 
w-plane which leads to a residue ("resonant contribution") when integration contours are shifted to 
the positive imaginary axis. The non-equilibrium polarizability gives only a resonant contribution 

a™(x', x) = i(e\iPf(x)^(x)^(x')^(x')\e) = J ^ e -M*'-*) 2 m 6(co - to 11 ) . (21) 

To illustrate the link to the retarded polarizability, we form the combination 

a e M, x) - a e M, x) = - [ —e^V-*) ( 1 + 1 

11V ' ' m ' ; J 2vr \oj n + u + iO+ uj n -uj-i0+ 

^e^K^r,), (22) 

which has poles in the lower half-plane only, as it should. For ground-state atoms, replace oj n o —lo u 
in these expressions, cf. Eq. Q for af™. 
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B Green functions of the free electromagnetic field 



We use the following Keldysh-Green functions for the photon field. To illustrate the formalism, we 
allow in this appendix for a thermal state with inverse temperature /3 whence 

D ll (x',x)=i(T{E(x')E(x)}) /3 = J ^e— (*'-')^ 11 (o;,r / ,r) , 

D 12 (x',x)=i(E(x)E(x'))p = i J (''- t )2n( W )Im[^(a;,r',r)] , 

D 21 (x',x) = i(E(x')E(x)) p =-i J C'- t )2n(-a;)Im[ J D i? ( W ,r / ,r)] . 

Note that the thermal occupation number becomes n(±u;) — >• ^fdfau) at zero temperature. The 
Fourier transforms D\2{uj) and L>2i(w) are then only supported by negative / positive frequencies, 
respectively. In the rest of the paper we always use this limit. The general case given above follows 



from the fluctuation-dissipation theorem [19 30 . The link to the retarded and advanced Green 



functions is provided by the relations Dr = Dx\—D\2 and Da = E)\\—D%\. The Feynman propagator 
in the frequency domain is, therefore, 

-Dii(w, r', r) = cwwo = Re[D^(oj, r', r] + i coth (/3w/2) Im[Z)^(w, r', r)] . (23) 

The retarded Green tensor in a linear and isotropic medium depends only on the difference r' — r 
and is given by 



Ayr /j 2 Akr 

D r (lo,t,0)= — 5(r)t + 



-,2 



where r = r/|r|, and k = y / s(u)){i(u>)uj/c (Im[k] > 0) is the wave vector in the medium. Obviously, 
[Dr(uj, r, 0)] 2 oscillates at half the medium wavelength, while \Dr(uj, r, 0)| 2 does not. In the paper, 
Dr(uj,tb,ta) is evaluated by contracting the tensor above from left and right with the transition 
dipoles d B and d A . 
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